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ABSTRACT 


Due to increasing concerns about CO 2 emissions and the shortage of fossil fuels, 
renewable energy has become a major topic in economic discussions. One renewable 
source is energy that can be extracted from the wind. This thesis covers the basics of 
using a doubly-fed induction generator (DFIG) to convert the mechanical energy of the 
wind into useful electrical power that can be used to supply electricity to any grid. 
Implementation and simulation results are analyzed in this research. The design 
implements digital four quadrant control of a DFIG with a direct current (DC) machine 
serving as the prime mover. Digital control of voltage, current, and frequency in the rotor 
windings is accomplished using a Voltage Source Inverter (VSI), while the stator voltage 
and frequency is maintained by the grid. Simulation is accomplished using Matlab and 
Simulink software. The simulations are verified with lab hardware. 
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EXECUTIVE SUMMARY 


Due to increasing concerns about CO 2 emissions and the finite supply of fossil 
fuels, renewable energy has become a major topic across the globe. Wind energy 
generation has attracted much interest in the last few years. Large wind farms have been 
planned and installed in various locations around the world. Many of these wind farms 
are based on the Doubly Fed Induction Generator (DFIG) technology with converter 
ratings around 30 percent of the generator ratings [1]. DFIGs have some advantages over 
synchronous and induction generators when used in wind farms, such as variable speed 
operation, active and reactive power control, and lower converter cost [1]. 

The objective of this research is to simulate and implement a wind energy 
conversion system via a DFIG. The control strategy utilizes a Field Programmable Gate 
Array (FPGA) and a Voltage Source Inverter (VSI) to implement digital four quadrant 
control of the rotor windings of a DFIG. Matlab/Simulink is the software resource used 
to design, simulate, and predict the system behavior before actual implementation. A 
Direct Current (DC) machine is used as mechanical load representing the wind energy 
source during the implementation. 


DFIG 



Figure 1. Simplified Circuit of DFIG Control for Wind Energy Conversion. 
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As seen in Figure 1, the voltage source inverter converts a 3-phase AC source to a 
controlled DC bus voltage. The DC voltage is then inverted into a 3-phase waveform via 
a digitally controlled VSI. The digital control for the Insulated Gate Bipolar Transistor 
(IGBT) network is implemented via a FPGA. The FPGA interface samples the voltage 
and current of a DFIG and speed of the DC drive through an analog to digital (A/D) 
converter. An algorithm then calculates the phase and amplitude of the voltage and 
current needed at the terminals of the rotor circuit. By controlling the phase and 
amplitude of the currents in the rotor windings various modes of operation and control of 
the DFIG can be accomplished. Real power flow and reactive power flow can be 
controlled, as shown in Figure 2, where the variables Vs and h represent the stator voltage 
and stator current. 


Unity Power Factor 



36 Degree Power Factor 



Figure 2. Power Factor Control of the stator circuit is achieved by changing the rotor 

currents. 
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A DC machine is used as the prime mover to simulate the wind. By manually controlling 
the torque produced by the DC machine, variable wind load conditions can be 
represented. As input torque is increased, the speed of the rotor changes and the position 
is fed into the algorithm of the FPGA. This corresponding position is then used to define 
the current in the rotor windings. The increase in input torque increases the output power 
produced by the system. See variable torque operation in Figure 3 where the stator 
current increases when the input torque increases. The variables Vs, h, and Ps represent 
the stator voltage, stator current, and average instantaneous output power for one phase 
multiplied by 3 from the stator. 


Step Change in Torque (1 N.m) 




200 



-0.01 0 0.01 0.02 0.03 0.04 0.05 


Figure 3. Step Change in Torque after steady state is reached (stator current, stator voltage 

and average instantaneous stator power). 
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The rotor current can also be used to control output power from the generator during 
steady state or transient conditions. An increase in the rotor current during steady state 
condition will cause a change in the stator current, which results in a change in the output 
power. Figure 4 shows a step change in rotor current and how it affects the output power 
for a given load on the shaft. The output power has slightly decreased because the 
increase in rotor current has caused the power factor to decrease. The variables h and Ps 
represent the rotor current and average instantaneous output power of one phase 
multiplied by 3 from the stator. 


Step Change 1-2 Amps 




Figure 4. Output Power vs. Step Change Rotor Current after steady state is reached. 

The goal of this research was to design an electrical interface for a wind energy 
conversion system. This thesis covers the simulation and hardware construction of a wind 
energy conversion system. The simulation results will be used to validate the actual 
implementation. Future research will investigate control schemes to improve on 

efficiency, effectiveness, and possibly grid fault analysis. 
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I. INTRODUCTION 


A. PURPOSE 

The mission at the Naval Postgraduate School (NPS) is to provide students with 
the highest quality and most defense-relevant graduate education available in computer 
and electrical engineering. The United States President’s Agenda on Energy and the 
Environment encourages the development of energy sources that will increase national 
security by decreasing our dependence on foreign oil, ensure that 10 percent of our 
electricity comes from renewable sources by 2012, and 25 percent by 2025 [2]. Wind 
energy is a very promising renewable energy source of the future and is part of the 
president’s overall strategy to develop a diverse portfolio of domestic energy supplies. 
This thesis supports the NPS mission and one topic discussed in the President’s Agenda 
on Energy and the Environment by using a wind energy conversion system as a means of 
increasing national security by harnessing one of nature’s renewable energy sources. 

B. RESEARCH OBJECTIVE 

The goal of this thesis was to build a computer simulation of a Wind Energy 
Conversion System, based on mathematical equations, through the use of Mathworks’ 
Simulink software. Once the simulation is complete, a fully functioning model of a Wind 
Energy Conversion System is built, so that various tests on the simulation and the 
physical system can be performed to compare the simulation results with the physical 
model. Based on the successful results of the system, we lay the groundwork for future 
research in the NPS power labs in the field of renewable energy. 

C. APPROACH 

The first step is to generate a computer simulation of the system, shown in Eigure 
1. There are numerous ways of simulating a DEIG, but the one used for this thesis is 
covered in [3]. In particular, a computer model of a DEIG based on mathematical 
equations is simulated in the rotor reference frame. A current control algorithm is 

designed in Simulink along with the Xilinx toolbox that controls the rotor currents to 
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produce the required voltage and frequency needed on the power grid. The wind is 
simulated via a torque input block. This model allows for representation of various modes 
of operation, which will be useful for future research, but for this thesis we concentrate 
on variable load conditions on the rotor, which simulate changes in the wind speed. With 
a working simulation and the ability to control rotor currents, various modes of operation 
and control are analyzed through the Simulink software. 

The physical system consists of a DFIG, VSI, Student Design Center (SDC) [4], 
encoder, and a DC drive used to simulate changes in wind speed. The DFIG is connected 
to the grid via the stator windings and to the VSI via the rotor windings. The SDC is the 
mechanism by which the control algorithm is implemented. The SDC measures the 
applicable system parameters and outputs the appropriate control signals to the VSI to 
control the rotor currents in a fashion that allows for a specific mode of operation. In this 
research, variable wind speed control, and active and reactive power control are 
demonstrated. 

D. THESIS ORGANIZATION 

This thesis is organized as follows. Chapter II covers the theory of induction 
machines. The equations that represent voltage, flux linkage, and reference frame 
transformations are discussed. The testing procedure used to measure and calculate the 
machine parameters used for simulation are also discussed. Lastly, Chapter II explains 
how the voltage and flux equations and machine parameters are used to model the 
system. Chapter III is dedicated to the control algorithm. It explains how an algorithm is 
formulated based on the equations derived in Chapter II and are manipulated to control 
the simulation and the physical system. Chapter IV displays the results of the working 
system and compares the actual results to the simulation results. Chapter V discusses the 
conclusions and the continuation research objectives for this thesis. Appendix A, B, C 
and D contain Matlab files, data sheets for equipment used, a full schematic layout of all 
Simulink systems, and transformation derivations used in this thesis. 
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II. DFIG SIMULATION 


A. MATHEMATICAL REPRESENTATION OF DFIG 

In order to simulate the Wind Energy Conversion System a proven method is 
needed to represent the characteristies of a DFIG. The equations used to charaeterize the 
DFIG used in this researeh are obtained from [3]. The rest of this seetion is devoted to the 
mathematical equations that describe a DFIG. 

The voltage equations to represent a 3-phase induction machine can be expressed 
as follows, 

^abcs ^ ^abcs P '^abcs ^ ^ 

^ abcr ^abcr P ^abcr 

where v, r, i, and X respectively refer to the voltage, resistance, current, and flux linkage 
of the phase windings. The subscripts a, b, and c refer to their phase component. The 
subscripts r and 5' refer to the stator and rotor windings. The term p represents the 

derivative ^ . The flux linkages shown in (1) for a linear magnetically coupled circuit 
are expressed as 

^abcs ^ 

_^abcr _ ^sr 

where L refers to the inductance in the respective winding. The subscript sr represents the 
mutual coupling between the stator and rotor windings. Finally, to complete equations 
that describe an induction machine the winding inductances are defined as 

T +T -if -if 

Hs ^ ‘-'ms 2 ‘-'ms 2 ‘-'ms 

— if T a-T —if 

-if -if f +f 

2 ‘-'ms 2 ‘-'ms ‘-is ^ 
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2 ^mr 


-ih 
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2 ‘^mr 
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cos 0^+^ 

cos 9^-^ 

Or-^-T 

cos 6^ 
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Or+^-f 

0r+^ 

cos 6^-^ 

cos 6*,. 


( 5 ) 


where the variables Lis and Lms are the leakage and magnetizing inductances of the stator 
windings, Lir and Lmr are for the rotor windings. The variable Lsr is the amplitude of the 
mutual inductance between the stator and rotor windings and 0r is the angular position of 
the rotor. 

As a standard practice and for convenience, the rotor variables are referred to the 
stator windings by the appropriate turn ratio N. 






^abcr ^abcr 






( 6 ) 


( 7 ) 


( 8 ) 


( 9 ) 


The magnetizing and mutual inductances are associated with the same magnetic flux path 
and related by the following equations. 



( 10 ) 
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L' = 


K^rJ 


( 11 ) 
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Ir 




( 12 ) 


By defining L;=^4^and 1^ = 
referred to the stator windings as 




the rotor inductance variables can be 


LL = L 


cos 6^ 

cos 0^-^ 
cos 0^+^ 


cos 0^+^ cos 0^-^ 


COS0^ 


cos 0,-^ 


cos 0^+^ 

COS0^ 


(13) 
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a:+4 
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-1 T 
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IT 

2 ^ms 


IT 
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(14) 


Substituting (13 and 14) into (1 and 2), the flux linkage and voltage variables referred to 
the stator windings, can now be expressed as 


^abcs 




^abcs 


s 

sr 

2 ' 

_^abcr _ 


-1 

bn 

_1 

_}abcr _ 


(15) 


Vabcs = r, 

^abcs P ^abcs 

V'abcr = K 

^abcr P ^abcr 


(16) 


The voltages and flux linkages in (15 and 16) describe the behavior of a DFIG 
and have components that vary with time and rotor position. A change of variables is 
used to reduce the complexity of these differential equations. A change of variables that 
formulates a transformation of balanced 3-phase equations to a new reference frame may 
be expressed as 
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( 17 ) 


[/;«.]=c /*, 
[/;v]=^(; /.; 


abcr 


where the matrix for stator variables and for rotor variables is expressed as 



cos6' 

sin6> 

2 


cos 6-^- 
sin 6-^ 

2 


cos 6 + ^- 
sin 0 + ^ 

2 


(18) 



cos P 
sin p 

. 2 


cos p-^- 
sin P-^ 

2 


cos P + ^ 
sin + ^ 

2 


(19) 


The variable 6 in (18 and 19) is the position of the reference frame for the variables 
andy^ = 6>-6(.. The x in (17) represents the particular frame of reference that the 
variables are transformed while / represents the voltage, current, or flux linkages in the 
stator or rotor windings. In the remaining equations in this chapter the absence of a 
variable in the place holder of x represents an arbitrary frame of reference. 


= 


cos 6y-6^ 
-sin 6^-6, 


sin e^-e^ 0 

cos Oy-0^ 0 

0 1 


( 20 ) 


To transform between reference frames use (20) where y represents the new frame of 
reference. Utilizing a little trigonometry, and performing the transformation in (17) on 
each component in the matrices on the right hand side of (16), the voltage equations can 

now be expressed in the arbitrary reference frame as 
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^ds=^sL-^^ds+P^ds 

^Os=^shs+P^r 

v' = r'i' + (Z) — (Z) A' + pA' 

qr s qs r qs qs 

4 =< 4 - (0-0)^ ^'ds+PKs 

^Or ~ ^shs P^r (21) 

where co^ is the rotational speed of the rotor and (O is the speed of the reference frame 

that the variables are transformed in radians per second. The next step is to express the 
inductance component of the flux linkages in (21) as a reactance instead of as an 
inductance. This is convenient because in the next section (Induction Motor Test) the 
actual machine parameters are measured and calculated in ohms. The inductance is 

converted to reactance at a base frequency co^ of 377 radians per second and (21) is 
rewritten as 


0) p 

0) p 

^ds = ^Ads - W,s + — Wds 

0 ), 0 ), 

Vos=rsios+ — ¥os 


Kr = '■/C + 


co-co^ 


COu 


-¥dr 


p 
CO. 


qr 


''dr 


= ^hr-- 


co-co^ 


COi. 




qr 


0), 


Vo, = rj. 


Or +~¥L 

OJu 


( 22 ) 


The flux linkage component of (21) now has units of flux linkage per second and 
represented by the symbol y /. For ease of understanding the flux linkages in (21) and flux 
linkage per second in (22) are compared 
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^ — L] i L (i i ^ 

^ds ^Is^ds ^ms ^ds ^dr ) 

■^S ~ ^s^Os 

^qr ^idqr ^ms ^qs ^qr 
^dr ^ddr ^ms ^ds ^dr 
■^r ~ ^idor 

Wqs = ^Jqs + ^msdqs + ^qr) 

Wds = '^Jds + ^msdds +C) 

Qs ^Is^Os 

¥'qr = ^'idqr + ^ ms ^qs + C 
W'dr = ^'irhr + ^ ms hs + Cr 

V^Or=X;Ar- 


( 23 ) 


(24) 


Equations (21 and 22) are derived in this section and mathematically represent the 
behavior of a DFIG. The next section discusses how to measure and calculate the 
parameters of the actual machine used in this thesis. 
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B. INDUCTION MOTOR TEST 


In order to efficiently model the system in Matlab and Simulink, the physical 
machines parameters are used for the simulation. The DFIG, model 8231, used for this 
research is part of the Lab-Volt Electro-Mechanical System (EMS) and consists of a 4- 
pole, 0.25 horsepower, 120 volt, 60 hertz, wave wound rotor machine. The stator consists 
of a set of 3-phase wye-connected windings with turns Ns and resistance Rs. The rotor 
consists of a set of 3-phase wye-connected windings with turns and resistance R^. 

The induction motor test performed is carried out as in [5]. 

The steady state operating characteristics of the DFIG was analyzed using a per- 
phase equivalent circuit shown in Figure 5. The DFIG was converted to an Induction 
machine for the parameter test by short circuiting the rotor winding terminals in a wye- 
configuration. The variables , R^ and represent the voltage, resistance, and leakage 

reactance of the stator. The variables V ', R '^, and represent the voltage, resistance and 

leakage reactance of the rotor referred to the stator. Xms represents the stator magnetizing 
reactance and s is the slip. 


Rs 


Xls 


Rr' 


XW 



Rl 


l-s 


Figure 5. Induction Motor Equivalent Circuit (From [5]). 


1. DC Test 

The first parameter test ran was to calculate the stator winding resistance for each 
phase. This test consisted of applying a known DC voltage to the stator terminals and 
measuring the current flow through the windings. The DC test on the stator windings is 
useful because it does not induce a voltage in the stator or rotor windings and eliminates 
the need to measure or calculate any reactance caused by induced voltages. The only 
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circuit parameter limiting current flow is . The measured voltage and current values are 

inserted into (25) and the stator resistance Rs is calculated. This test is performed on each 
phase of the stator windings and the three values are averaged together. 

(25) 

^dc 


2. No Load Test 

A no-load test is performed to measure the rotational losses and other equivalent 
circuit parameters needed in the blocked rotor test. In this test, rated voltage and 
frequency is applied to the stator while the machine is running at no load. The input 
power, phase voltage, and phase current are measured for each phase and averaged. At no 
load, the slip is very small and the rotor branch impedance is very large. Therefore, all 
the real power is dissipated in the stator resistance. The machine at no load is very 
inductive, so the input impedance is expressed as 

|4,«l| (26) 

From this equation, we can solve for the quantity , which is used in the blocked 

rotor test. 


3. Blocked Rotor Test 

The blocked rotor test is performed to calculate the remaining circuit parameters. 
In this test, the rotor of the induction machine is blocked and a reduced voltage is applied 
to the terminals so that rated current flows through the stator windings of the machine. 
The input power, voltage, and current are measured and averaged. When the rotor is 
blocked, the slip is equal to one. As a result, the rotor branch impedance is very small and 
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the magnetizing impedance is ignored. Real power is only consumed by the stator and 
rotor resistances. The sum of the stator and rotor resistances can be solved as 


R,+R' = 




3/, 


|2 • 


The input impedance is expressed as 


Ik,* 


7 R.+K '+ ' 


(27) 


(28) 


where it is assumed thatX;^ =^ir- Substituting the result of (27) into (28) allows us to 
solve for the numerical quantities of and . From this point X^ can be calculated 

by substituting the stator leakage reactance value into equation (26). Table 1 below is a 
summary of parameters calculated from the DFIG used for this thesis. 


Parameters (0) 

Stator Resistance 

12 

Stator Leakage Reactance 

9.1 

Rotor Resistance r/ 

15 

Rotor Reactance X^ 

9.1 

Magnetizing Reactance 

126 


Table 1. DFIG Measured Values at 60 Hz. 

C. DFIG COMPUTER MODEL 

The equations (22 and 24) are the basis for simulating the DFIG. A block diagram 
of the simulation setup is shown in Figure 6. The model simulates the DFIG in the rotor 
reference frame. Although the equations that represent a DFIG have been defined, they 
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still need to put in a form that is useful for computer simulation. We also need to define 
the torque, power, and rotor speed equations, which are needed for simulating the 
machine. 



Figure 6. DFIG Computer Model in Block Form (After [3]). 

The first step is to solve the flux linkage per second equations of (22) for current. 
We do this by defining four new variables in terms of quantities calculated in the 
induction motor test 




ms qs qr 


aq 


Kj 


¥md=^ms L+i'dr =X. 


ad 




(29) 


(30) 
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where 


Y — Y 

^aq ~ ^ad 


" 1 1 1 ^ 


■+-+- 


Kj 


The currents are now defined as 


Z,. 
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(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


Now that the current equations have been defined, the flux linkage per second equations 
in (22) are rewritten as integral equations defined in term s of the four new variables from 
(29-31). The final form of the flux linkage equations that are used to simulate the DFIG 
shown in Figure 6 are represented as 


y^,s=(^b 


i 


CD K 

^qs - V^ds+YT 


COu 


! mq r qs 

^Is 


(38) 


Wds=0^b\ 


CO r 

^ds+ — W,s+— Wmd-Wds 

A,, 


OOu 


(39) 


13 



(40) 


=(Ob\ 

J 




X,. 


¥[r = (Ob \ 



f \ 

f v' - 

co-co^ 

^qr 

\ (Ob ^ 

J 


¥dr+^ Wm,-¥c 
^Ir 




f 

\ 

( 

f 

co-co^ 


(Ob\ 

^dr + 



J 


1 (Ob 

J 


Wdr+i;^ Wnui-Wdr 

^Ir 


V^r 


, f r r' , 

ar=(Obi 


(41) 
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(43) 


The rotor position, which is needed to transform the equations to the respective 
frame of reference, is 


0 ). 


f T-Tr P 
= \— - -—dt 

J 9 / 


(44) 


where and 7), represent the electrical torque generated in the machine and input torque 

on the shaft of the rotor. The variables P and J are the number of machine poles and the 
inertia of rotor. Once the speed of the rotor is obtained, the position is calculated as 


6^ = jo)^ dt. 

The electrical torque and the power equations are defined as 


(45) 
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D. SIMULINK IMPLEMENTATION OE DEIG 


The DFIG from Figure 6 is implemented with Simulink Software. The following 
Figures 7 and 8 show how (38-46) are implemented using the Simulink software to 


represent the DFIG. 



Figure 7. DFIG Simulink Simulation (From [6]). 

Figure 7 shows the blocks where the flux linkage per second equations are implemented, 
and how they are used to compute the rotor and stator currents. The respective equation 
number is annotated on each block along with the associated input and output. Figure 8 is 
a drill down into the block that represents (38) the stator flux linkage per second on the q- 
axis. As seen in Figure 7, there are four inputs and one output, which also serve as an 
input to the same block. Using the basic mathematical operators in the Simulink library 
(38) is represented as shown in Figure 8. The remaining blocks for each equation that 
represent the DFIG are built the same way and can be found in Appendix C. 
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Figure 8. Equation 38. 

In conclusion, this chapter described the methods in which to simulate the DFIG 
used in this thesis. The next chapter will discuss the methodology used to control the 
Wind Energy Conversion System. 
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III. CONTROLLER 


A. CONTROLLER TOPOLOGY 

As stated earlier, one of the components of the wind energy conversion system is 
the rectifier and VSI. The VSI and rectifier used for this research are a combined package 
manufactured by SEMIKRON and called a SEMISTACK. The SEMISTACK consists of 
a three phase rectifier and three phase inverter. The inverters inside the SEMISTACK are 
made of SKM 50 GB 123D IGBTs that are controlled by SEMIKRON SKHI 22 gate 
drivers. The VSI is used to control the magnitude, frequency, and phase angle of the 
currents in the rotor. The rectifier produces a constant DC voltage to the VS Is. The 
control for each IGBT in the SEMISTACK is implemented through the SDC. The SDC 
utilizes a Xilinx Field Programmable Gate Array (FPGA) chip and a control board that 
has various inputs and outputs that interface with the SEMISTACK and DFIG. One of 
the inputs is an 8-channel, 12-bit A/D converter used to measure the voltage and current 
signals from the DFIG. Another input samples the rotor speed. The pre-programmed 
FPGA chip, based on the sampled inputs, outputs the desired gate signals to the control 
board. The control board has BNC (Bayonette Neil-Concelman) connecters that connect 
to the SEMISTACK’s gate drivers. The gate drive signals activate the IGBTs in a fashion 
that will produce the desired 3-phase current in the rotor windings. The SDC has many 
parts and functions, which make it useful as a design and analysis tool. The procedure for 
using the SDC, architecture of internal components, and the benefits of using a FPGA for 
digital control of power systems are discussed in [4]. 

Figure 9 is a Simulink block diagram of the Wind Energy Conversion System. 
The diagram shows where the controller samples each input and the corresponding 
outputs that go to the DFIG. The controller is physically located inside the SDC for the 
hardware design. It is important to note that the same controller logic used for simulation 
is used in the hardware design; therefore, the input/output blocks shown in the simulation 
block diagram symbolize actual connections to the SDC from the grid and the DFIG. 
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Figure 9. Simulation of DFIG Control for Wind Energy Conversion. 


Figure 10 is a drill down of the Rotor Controller Block. The same logic blocks 
used in the simulation to control the DFIG model are used in the hardware configuration 
to control the DFIG with one exception. An additional set of logic blocks is added to 
implement Space Vector Modulation (SVM) which is one of many techniques that can be 
implemented on a VSI to create sinusoidal 3-phase waveforms and will be discussed 
more in Section B.7. For now, it is important to know that inside the rotor controller there 
is a separate SVM block that is not needed for simulation. Figure 10 shows how the rotor 
controller is used for the simulation and hardware design. The next section explains how 
each of the components in Figure 10 is used to control the rotor current. 


18 

























































Figure 10. Rotor Controller Block Diagram. 


B. CONTROLLER COMPONENTS 
1. Voltage Source Simulation 

The simulation required a voltage block to simulate the voltages on the grid. The 
design of the source is shown in Figure 11. The voltage block models the rated 
parameters of the DFIG used for this thesis. The voltage block produces three sinusoidal 

waveforms at 60 hertz with peak amplitude of 120*^/2and displaced equally 120 
degrees apart. The controller setup requires only two voltage inputs so the line to line 
values are extracted as shown. 



Figure 11. Voltage Source Simulation. 
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2. 


Transformation of Stator Voltage 


As seen in Figure 10, the stator voltages are fed into the rotor controller where 
they are transformed to the rotor reference frame. The transformation fabcs 

in (17) is slightly modified to transform the line-to-line voltages and to the rotor 

reference frame (See Appendix D). Figure 12 demonstrates how the transformation is 
performed with the Xilinx library set. The constant 85 in Figure 12 corresponds to a 30 
degree phase shift from the transformation in Appendix D. The sin and cosine blocks 
from the Xilinx library use a 10 bit lockup table to represent angles from 0-360 degrees. 
The angle 9^^ is obtained by taking the angle from arctangent of the transformed voltages 

and represents the electrical position of the currents in the rotor. This angle is used to 
transform the rotor currents to the synchronously rotating reference frame. The voltage d- 
axis component is inverted before taking the arctangent to account for the d-axis being 
180 degrees out of phase from the y-axis. 



Figure 12. Stator Voltage Transformation to Rotor Reference Frame. 

3. Transformation of Rotor Currents 

The rotor currents 4^^, shown in Figure 13, are also fed into the controller where 
they are transformed using (17) to the variables . This notation represents a 
transformation of the rotor currents to the synchronously rotating reference frame. This 
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transformation produces a constant value, if the currents are in steady state condition, 
which is useful in the next block. The constant 341.33 represents a 120 degree phase shift 
used in (17-19). 


341.33203125 






b^Cjb) 


b-'Cjb) 


b^Cjb) 


¥ 


: (bb) -I 


Figure 13. Rotor Current Transformations. 

4. Proportional /Integral (PI) Gain 

The transformation of balanced 3-phase variables to the qdOaxis will only 
produce components on the qd axis. Additionally, if the transformation is to the 
synchronously rotating reference frame, the q component will have constant amplitude 
equal to the amplitude of the phase variables and the d component will have amplitude of 
zero. This property is useful in the PI control algorithm for driving the steady state error 
to zero. Figure 14 contains the logic used for a proportional and integral gain controller 
designed with the Xilinx library block set. By transforming the rotor currents to the 
synchronously rotating reference frame two constants are produced as discussed above. 
Setting the reference iq to the desired output value and the reference id to zero allows for 
full control of the output. The Integral part of the PI controller is built with an 
accumulator so a reset is added to set the accumulator to zero when desired. The currents 
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in the rotor windings are controlled by adjusting the voltage at the rotor terminals; 
therefore, the output from the current PI controller is a voltage. 




Figure 14. PI Current Controller. 

5. Transformation of Rotor Voltages to Rotor Reference Frame 

£ 

The rotor voltages v® and are transformed from the synchronously rotating 

reference frame to the rotor reference frame as shown in Figure 15. This transformation 
utilizes equation (20), which is useful for transforming variables between frames of 
reference. The input to the VSI in Figure 10 is normalized to the value of the DC voltage 
produced by the rectifier on the SEMISTACK. This normalization takes place in the last 
logic block of Figure 15. 
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Figure 15. Transformation Between Reference Frames. 

6. Rotor Encoder 

The encoder used for this thesis is an MES20 (Type C) incremental shaft encoder. 
The encoder provides a digital input to the SDC. The control circuitry inside the SDC 
interprets the digital signals from the encoder and produces both rotor speed and position. 
The encoder software used in this thesis is discussed in detail in [7]. One modification, 
shown in Figure 16, was added to the design in [7], which was used to adjust rotor 
position and speed for a 4-pole machine. See Appendix C for complete encoder layout. 


From Encoder 





Figure 16. Rotor Encoder Addition. 
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7. Space Vector Modulation 

One well-known way of creating 3-phase waveforms from a VSI is via SVM. 
SVM is a form of Pulse Width Modulation (PWM) in which an algorithm involving 
space vectors is used to control the on and off times of pulsed signals. The pulsed signals 
then drive the input signals to a VSI allowing the user to create any magnitude and 
frequency of signal desired. The SVM approach utilized in this controller is obtained 
from [7]. 

The VSI produces a three phase output voltage that is controlled by turning on 
and off the six IGBTs on the VSI diagram in Figure 17. The on and off switching states 
are based on the SVM Hexagon also shown in Figure 17. The qd axis is overlaid on the 

SVM Hexagon and represents the axis of the variables and , which are shown 

entering the SVM block in Figure 10. The magnitude of the vectors and is equal to 

the vector V * and 6 is equal to the arctangent of the two vectors. 



Figure 17. VSI and SVM Hexagon (From [7]). 


There are six sectors on the hexagon and eight possible on and off states that are used to 
produce the vectorsV*. The vectors Vj and when summed are also equal to the 

vectorV*. In Sector I, VJ andV 2 » correspond to the states (p,n,n)and(p,/7,n). The 
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pand n correspond to the positive or negative IGBT bus signals for the respective 
phases andV^. The zero vectors, which are represented by the states {p,p,p) 

and(n,n,n) are not shown but represented as vectors into and out of the page. The 
magnitude of the vectors and V 2 correspond to the amount of time spent on the 
switching states in each sector. The duty cycles 7^ and are the times spent on each 
cycle for the vectors andV^ ■ The total time spent for one switching period is 7^. 




^2 = 


2Ty^ 


dc 


37; 

3r 


(48) 


(49) 


Applying the law of sines to any one of the sectors in the SVM Hexagon will produce 


2y* 


3) 


y, 


^/3 sin 60°-^ sin 0 


(50) 


Substituting (48 and 49) into (50) is used to find the time of the duty cycles for the 
vectors andy 2 . 


V*J3 

60'’-0 

i X 7 j 


y 


dc 


T ■ a 

u =-sin 9 


dc 


T^=T,+T2+T, 


(51) 


(52) 

(53) 


SVM has the ability to minimize switching loss and harmonic distortion by controlling 
the order in which the states are applied [8]. The switching pattern for each state and 
sector used in this thesis is shown in Table 2. The SVM algorithm is implemented inside 
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the SDC center where it is oversampled to converge on the performance of a true 
continuous sinusoidal signal. The digital implementation of the SVM scheme used for 
this thesis is discussed in detail in [9] and the schematic diagrams can be found in 
Appendix C. Figure 18 is a block diagram of the functions for each process in the SVM 
block. 
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Table 2. Space Vector Modulation Switching Pattern (From [9]). 
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Figure 18. SVM Digital Implementation (From [9]). 

In Conclusion, this chapter described the methodology used to control the DFIG 
for simulation and hardware design. The next chapter discusses the test performed to 
verify the operation of the system. 
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IV. RESULTS AND ANALYSIS 


A. VOLTAGE AND CURRENT COMPARISON 

Several tests were performed to verify the operation of the wind energy 
eonversion system. The first test is used to verify that the system ean produee rated stator 
voltage and frequeney on the grid. Figure 19 is a eomparison of the physieal system to 
the simulation with the rotor current set at lA. The simulation closely matches the 
physical system and produces a 170V peak to peak signal at 60 Hz. The variables Vs, and 
Vs,Sim represent the stator voltage and stator simulation voltage. The variable L represents 
the rotor current. The variable Pout represents the average instantaneous output power 
from one phase multiplied by 3 delivered from the stator terminals to the grid. This 
notation is followed for all the graphs in the results section. 


Comparison of Stator Voltages 



Figure 19. Comparison of Grid Stator Voltage to Simulation Source. 


29 

























































The second test was to verify the operation of the system for various changes in 
rotor current. The control algorithm allows the user to specify the amplitude of the rotor 
current. The benefits of controlling the amplitude of the current are discussed in the next 
section. This test is just to prove that the system functions similarly to the simulation for 
changes in rotor current. The results are shown in Figure 20. 

Generator Stator Currents 



0 

T3 

Q. 

E 

< 



Time (s) 

Figure 20. Comparison of Stator Currents to the Stator Simulation Currents. 

From Figure 20, one can see that as the rotor current is increased in amplitude, the stator 
current reduces in amplitude. The simulation does not produce as much change in stator 
current amplitude as the physical system for equal changes in rotor current. This is not a 
problem; the system still produces voltage and current at 60Hz. 
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B. POWER FACTOR AND OUTPUT POWER COMPARISONS 


The next test demonstrates the usefulness of having the ability to control the 
amplitude of the rotor currents. A Wind Energy Conversion System can only output as 
much power as the wind inputs into to the turbine blades. A specific amount of torque on 
the turbine blades will generate a certain amount of apparent power in the machine. 
Controlling the amplitude of the rotor currents allows the user to determine how much 
real power is abstracted from the total amount of apparent power generated in the 
machine. This technique is known as power factor control and shown in Figure 21. The 
voltage and current values below are normalized to the DFIG ratings. 


Power Factor Comparison 



Power Factor Comparison 



Power Factor Simulation Comparison 



Power Factor Simulation Comparison 



Figure 21. Power Factor Comparison. 

Figure 21 compares the power factor adjustment for the physical system with the 
simulation. The rotor current is adjusted from 1 to 3 amps and the stator voltage and 
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currents are graphed. One can see that the angle between the voltage and current are 
almost in phase, as the rotor current is shifted from 1 to 3 amps. 

Figure 22 is a display of how the power factor affects real power delivered to the 
grid for changes in rotor current. The power factor for rotor currents of lA, 2A, and 3A 
respectively are calculated as 0.3, 0.71, and 0 .94 for the physical system and 0.28, 0.69, 
and 0.97 for the simulation. 
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Figure 22. Comparison of Output Power for Changes in Rotor Current 

The last test is used to verify that the system can operate under various changes in 
load. A DC machine is used to simulate a change in wind speed by increasing the amount 
of torque delivered to the DFIG shaft. A step change of 1 Newton is applied by the DC 
drive, which is equivalent to 175 watts at rated speed. The results, shown in Figure 23, 
prove that the system is capable of delivering constant power to the grid as wind speed 
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changes. The physical system produces slightly less output power than the simulation due 
to mechanical coupling losses from the DC drive to the DFIG. 
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Figure 23. Comparison for Step Change in Torque 

This chapter covered the results for each test performed that compared the 
operation of the simulation to physical system. The next chapter will discuss the results 
and conclusions. 
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V. CONCLUSIONS AND SUGGESTIONS 


A. CONCLUSIONS 

A Wind Energy Conversion System was designed, simulated, and constructed 
using the methods discussed. A series of tests were conducted to validate the design. The 
system was designed to work under various load conditions, which in this scenario 
simulates changes in wind speed. The ability to control the rotor currents provides the 
user with functionality to adjust power factor thereby increasing the efficiency of power 
delivered to the grid. The response of the Wind Energy Conversion system was compared 
to the simulation. The results show that the physical system behaved as predicted by the 
simulation. 

A successful Wind Energy Conversion system demonstrates the ability to 
transform mechanical energy delivered by the wind into electrical energy that can be used 
to power any electrical grid. The Wind Energy Conversion System successfully 
transformed mechanical torque on the shaft into electrical power. 

B. FUTURE RESEARCH OBJECTIVES 

This thesis dealt with simulation and operation of a Wind Energy Conversion 
System. Future research will be on various control strategies to improve overall 
efficiency. 

One area of research currently being explored is the ability to extract power from 
both the rotor and stator windings when the machine is operating at super-synchronous 
speed. This design will require two VSI and control software that will allow for bi¬ 
directional flow of current through the rotor windings for sub-synchronous and super- 
synchronous operation. Delivering power from the rotor and stator windings will improve 
the overall efficiency of the system. 

The other area of research currently being explored is a continuation of 
controlling the power factor angle. The rotor side converter will be used to control active 
and reactive power of the DFIG by controlling the amplitude of the rotor currents. The 
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amplitude of the rotor currents will follow a tracking characteristic that adjusts generator 
speed for optimal power generation depending on wind speed. 
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APPENDIX A: DATASHEETS 


SEMISTACK - IGBT 



SEMITRANS Stack'> 


Three-phase rectifier + 
inverter with brake 
chopper 

SEMITEACH - IGBT 
SKM 50 GB 1£3D 
SKD 51 
P3^50F 

Features 

* Muhti-functian IGBT converter 

* Transparent enclosure to allow 
visualization of every part 
IP2x protection to minimize 
safety hazards 

* External banana^BNC type 
connectors for all devices 

* Integrated drive unit offering 
short-circuit detection/cut-off, 
power supply failure detection, 
interlodk of IGBTs -i- galvanic 
isolation of the user 

* Forced-air cooled heatshk 

Typical Applications 

* Education: One stack can 
simulate almost all existing 
indusfrbl applications: 

- 3-phase inverter+brake chopper 

- Buck or boost converter 

- Single phase inverter 

- Single or 3-phase rectifier 

Phola non-conHadEual 


Circuit 

Inns (A) 

Vac ^ 

Types 

Beci 

30 

440/750 

SEMITEACH - IGBT 


Symbol 

Conditions 

Values 

Units 

Irma 

no overload 

30 

A 


IGBT - 4x SKM 50 GB 1230 



VCES 


1200 

V 

Vce[Sati 

Vgjes 

lc= 50A. Vge= 15V. Cfi|] level: Tj= 2£{125)'X; 

27(3.5} 

iZO 

V 

V 

Ic 

ICM 

T.^= 25 {BOfC 

T^= 25 {BOfC: tp= 1ms 

50(40} 

100(80} 

A 

A 

^iT^rran:) 

Hectilier - Ik SKD 5t/14 

without filter 

3x4B0 

V 

with filter 

3*380 

V 

CgcflH 

^rJCmBK 

DC Capacitor bank - Electrolytic 2i 2200pF/40CiV 

total equivalenE capacitance 

max. DC voltage applied to the capacitor bank 

1100/800 

750 

pF/V 

V 


Drrver-4x SKHI 22 



Power 

supply 

Current 


0/15 

V 

consump 

lion 

max: per driver 

16 

mA 

Ttiermal 

Irf] 

Normally Open type (NO) 

71 






_ n n n n n n n 

r*i' irtt 


nananana 


m 




Jf ft If r; 


D0OO 


-3D3- 

-3Dg- 
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Model 8211 DC Motor/Generator 



This machine can be run independently as a DC motor 
or a DC generator. The armature, shunt field, and series 
field windings are terminated separately on the faceplate 
to permit long and short shunt as well as cumulatively 
and differentially compomded mota and generator con¬ 
nections. This machirte is fitted with exposed movable 
brushes to allow students to study the effect of armature 
reactron and commutation while the machine is operatmg 
under load. An independent, circuit- breaker protected, 
sturft-field rheostat is mounted on the faceplate for motor 
speed control or generator output voltage adjustment. 


Model 8231 Three-Phase Wound-Rotor Induction 
Motor 



Each phase of the stator windings of this motor is inde¬ 
pendently terminated and identified on the faceplate to 
perniit operation in either delta or star (wye) configura¬ 
tion. The rotor windings are brought out to the faceplate 
via external slip rings and bmshes. This machine can be 
used as a wound-rotor induction motor, phase shifter, 
single-phase variable coupling transformer, three-phase 
transformer, selsyn cofTtrol, frequency converter or 
asynchrorwus induction generator. The speed of this 
machine can l3e controlled through the use of the Three- 
Phase Rheostat (Model 8731). 


SPECIFICATIONS 


Model B211 DC Motof/<3eiHf3lor 

120^09 V. 60 Hz 

22)07390 V- 50 Hz 

240W15 V- 50Hz 

Power Requrement 

12D.'20aV 

22Di330V 

246^15 V 

Rat ng Motor Output Power 

General or Output Power 

Armature Voltage 

Bhuil Retd Voltage 

Full Load 

Ful Load Motor Current 

Full Load Generator Current 

175W 

12flW 

HOW 

i2aw 

t 2i3 V- DC 

22aV- DC 

240 V- DC 

t20V- DC 

220V- DC 

240 V- DC 

1BOO nhiin 

1500 dmin 

1500 r7mii 

2BA 

1.3A 

1.1 A 

1 A 

D.5 A 

6.5 A 

Physical Charaderistics Dimensions (H Jt W Jt D) 

NIet Weight 

3Mx2ei K 440 rrm {12.1 3e 11.5 k 17.3 ii) 

14.1 kg (31 lb) 

Model 9221 Four-Pole Squirrel-Cage Induotiun Mator 

1207209 60 Hz 

2207390 V- 50 Hz 

240W15 V- 50He 

Power Requrement 

12D.'20aV 

220.330 V 

246415 V 

Ratng Output Power 

Stior Voltage 

Full Load 

Full ILoad Current 

175W 

12a'2Da V. 3rghase 

22D!'3B6 V. S-pbase 

240415 V. 3-phase 

1670 n^in 

1350 drrin 

1395 r7mri 

1.2 A 

0.52 A 

0.46 A 

Physical CharacSerislics Dimensions (H x W Jt D) 

NIet Weight 

309x2ei X 440 rrin { 12.1 x 1 1.5 x 17.3 in] 

13.5 kg (29.7 lb} 

Model 9231 Three-Phase Wound-Rolor Inductiofii Motor 

12Ba09V- 60 Hz 

2207390 V- 50 Hz 

240415 V- 50 Hi 

Power Requiernent 

12Di'2MV 

22D330V 

240415 V 

Ratng Output Power 

Sttor Voltage 

Rotor Voltage 

Full Load GpffiU 

Full ILoad Current 

175W 

m^2Da V. 3phasB 

22D.'3BD V. 3-pbase 

240415 V. S-phase 

6D!'104V. 3-phase 

11DM9D V. 3-pbaBe 

120)203 V. 3-Fhase 

1506 nrtnin 

1240 dmin 

1315 r7mri 

1.3 A 

0.53 A 

0.4BA 

Physical Charaderislics Dimensions (H Jt W Jt D) 

NIet Weight 

303 x 201 X 440 mn 02.1 x 11.5 x 17.3 in) 

14 kg (3D.a lb) 























































POWER SUPPLY 
MODEL 8821 


The Power Supply provides fixed and variable AC 
and DC voltage sources, all terminated by color-coded 
4 mm safety sockets. Independent circuit breakers, reset 
at the front panel, protect the input to and output ilfom the 
Power Supply. Indicator lamps monitor the presence of 
input voltage in each phase. When a phase leg of the 
site's power service is out, the lamp goes off to reflect this 
condition. 


A voltmeter, connected through a selector switch, 
nvonitors the variable AC and DC outputs and fixed DC 
output A 24 V AC output provides a low-voltage supply 
required to operate other EM S equipment such as meter¬ 
ing modules and modules used in the Power Electronics 
Training System. 


SPECIFICATIONS 


Model SBi21 - Power Supply 

120/2QaV - 60 Hz 

220raa0 V - 50 Hz 

244W415V-S0 Hz 

Inpul Line VDltage 

Line Cu^Ten^ 

Service Installatkm 

12D.'2{JBV 

22D!'3BD V 

240.410 V 

15 A 

IDA 

2D A. 3-phase, 5 wires, star i[Wye>rMnnerved, including neutral and 
ground 

Outputs Rjced AC 3-Phase 

Variable AC 3-PhasB 

Variable DC 

Fixed DC 

Lew Power AC 

115 A 

22{]L!3flOV- 10 A 

240.'415V- 1UA 

0-12a.20aV-5 A 

D-220.'3aO V - 3 A 

D-240.415V-3A 

0-120 V-a A 

■D-22D V - 0 A 

0-240 V-5 A 

120 V-2 A 

220 V - 1 A 

240 V- 1 A 

24 V - 3 A 

Wall Outlet included 

I^MAL21-20 

NEMA L22-20 

CLIPSAL 55SO02O 

Power Cord 

3 m,i:iDft:i 

Physical Characters tics Dimensons (H x W x D) 

Nel Weigh! 

303 X 287 X 500 mm (12.1 x 1 1.3 x 19.7 in) 

ia.4 kg i;4D.5lh;i 
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APPENDIX B: MATLAB M-FILES 


A. MATLAB INITIAL CONDITIONS FILE 


DblFedlnd.m 


clc,clear 
Vdc=150; 

DC_Bus_Prot=2 0 0 ; 
omega_b = 2’^pi’^60; 
oversample=2; 


%Bus voltage maintained for Capacitor on VSi 
%Secure firing circuit at 200Vdc. 

% Base Frequency in radians/second 
% Determine oversampling rate in SVM block 


pulsect = 1800/oversample; 


step_ct=l; 
clock_freq=2 5e6 ; 
tstep=.0004; 


only to speed up simulation. 


%tstep=step_ct/clock_freq; %When compiling software us this line 

%More constants used for simulation 
twopibyS = 2’^pi/3; 
poles = 4; 

polesby2 J=poles’^ 12/2 / . 08 9; 

TL= . 2 5’^7 4 6/ (2/poles’^omega_b) ; 

%Calculated constants from induction motor test. 



Xls =9.1; 
Xm =126; 
Xlr =9.1; 


%D= (Xls+Xm)(Xlr+Xm)-Xm"^2; Defined in book for per unit values of 
inertia. 

rsbyXls = rs/Xls; 
rrbyXlr = rr/Xlr; 

Xaq = 1/(1/Xm+l/Xls+l/Xlr); 

Xad = Xaq; 

XaqbyXls = Xaq/Xls; 

XaqbyXlr = Xaq/Xlr; 

XadbyXls = Xad/Xls; 

XadbyXlr = Xad/Xlr; 

XaqbyXm = Xaq/Xm; 

XadbyXm = Xad/Xm; 

psi_qsic=0; 
psi_dsic=0 ; 
psi_qric=0; 
psi_dric=0 ; 

omegar_ic =omega_b; %Set the initial condition speed of rotor to 1800 
rpm. 

Kp_i=.231; %current Proportional gain 
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Ki i=415.7^1.042; 


%Current Integral gain 


%Used for Mealy State Machine in A/D Current simulation. 

F_mat = [0 0 0 1;1 1 2 0;2 2 3 0;3 3 0 0]; 

0 mat = F mat; 


%Everything below this line is used for the encoder controls 
%fvtool(outputl,output2) 

t_square=2e-4/10; %2e-4 produces 1500RPM...4e-4 produces 750RPM 
%8e-4 = 374RPM...10e-4=300RPM 

%there is 400 1/2 pulses per rotation therefore: 

%60sec/(400^^ l/2’^t_square) = correct RPM 
RPM=60/ (200’^t_square) ; 

Hz=RPM/60; 

sw_freq=15000; 

% sw_counter=round(f_clock/sw_freq-mod(f_clock/sw_freq,10)); 

%Counter for sawtooth for switching modulo 10 used so step_ct can be 10 
Total_Rotations=6; 

%Ctr=[l:2"8]/2"8/tstep/400^50;%16, 000,000 
%reciprocal=l./Ctr;%l/16,000,000=244.1406 
Ctr=[l:2"ll]; 

reciproca1=360/30/1. /Ctr;%360/30/l./Ctr 
% [Pk, Nk, Ok] 

output_vec=[0; ... %Nothing is occurring with [0,0,0] 

5; . . . %3Motor is turning in the CCW direction and crossing 

the 

%zero pt with [0,1,1] 

0; ... %0Nothing is occurring with [0,0,0] 

4; ... %2Motor is turning in the CCW direction with [0,1,0] 

0; ... %0Nothing is occurring with [0,0,0] 

3; ... %5[1,0,1] indicates CW rotation and resets rotation 
%counter 

0; ... %0Nothing is occurring with [0,0,0] 

2; ... %4[1,0,0] indicates CW rotation 

3; ... %5[1,0,1] indicates CW rotation and resets rotation 
%counter 

0; ... %0Nothing is occurring with [0,0,0] 

2; ... %4[1,0,0] indicates CW rotation 
0; ... %0Nothing is occurring with [0,0,0] 

5; ... %3[0,1,1] indicates CCW direction and resets rotation 
%counter 

0; ... %0Nothing is occurring with [0,0,0] 

4;...%2 Motor is turning in the CCW direction with [0,1,0] 
0]; %0 Nothing is occurring with [0,0,0] 


next 

output=next; 


[ 0 , 0 , 1 , 0 ; 1 , 0 , 1 , 0 ]; 
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B 


MATLAB M-FILE USED FOR SPACE VECTOR MODUALTION 


cverflowS.m 

function [sectorl, sector2, sectorS, sector4, sectors, 
overflows(x) 

%gain = xfix({xlUnsigned,10,7},2.359296/3);%for 60 hz 
gain = xfix({xlUnsigned,10,7},2.359296); %for 180 hz 
% t emp v=g a i nX ; 
tempv=x; 
if tempv<=171-1 

sectorl=xfix({xlBoolean}, 1) ; 
sector2=xfix({xlBoolean},0); 
sector3=xfix({xlBoolean},0); 
sector4=xfix({xlBoolean} , 0) ; 
sector5=xfix({xlBoolean} , 0) ; 
sector6=xfix({xlBoolean},0); 
z=xfix({xlUnsigned,10,0},tempv); 
elseif tempv<=2171-1 

sectorl=xfix({xlBoolean}, 0) ; 
sector2=xfix({xlBoolean}, 1) ; 
sector3=xfix({xlBoolean},0); 
sector4=xfix({xlBoolean}, 0); 
sector5=xfix({xlBoolean} , 0) ; 
sector6=xfix({xlBoolean} ,0); 
z=xfix({xlUnsigned,10,0},tempv-171); 
elseif tempv<=3’^ 171-1 

sectorl=xfix({xlBoolean}, 0) ; 
sector2=xfix({xlBoolean} , 0) ; 
sector3=xfix({xlBoolean}, 1) ; 
sector4=xfix({xlBoolean},0); 
sector5=xfix({xlBoolean}, 0); 
sector6=xfix({xlBoolean} ,0); 
z=xf ix ({xlUnsigned, 10,0}, tempv-2171) ; 
elseif tempv<=4171-1 

sectorl=xfix({xlBoolean},0); 
sector2=xfix({xlBoolean}, 0); 
sector3=xfix({xlBoolean} , 0) ; 
sector4=xfix({xlBoolean} , 1) ; 
sector5=xfix({xlBoolean},0); 
sector6=xfix({xlBoolean},0); 
z=xf ix ({xlUnsigned, 10,0}, tempv-3’^ 171) ; 
elseif tempv<=5’^ 171-1 

sectorl=xfix({xlBoolean},0); 
sector2=xfix({xlBoolean} , 0) ; 
sector3=xfix({xlBoolean} , 0) ; 
sector4=xfix({xlBoolean} , 0) ; 
sector5=xfix({xlBoolean},1); 
sector6=xfix({xlBoolean} ,0); 
z=xf ix ({xlUnsigned, 10,0}, tempv-4171) ; 

else 


sector6. 


sectorl=xfix({xlBoolean} , 0) ; 
sector2=xfix({xlBoolean} , 0) ; 
sector3=xfix({xlBoolean},0); 
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end 


sector4=xfix({xlBoolean},0); 
sector5=xfix({xlBoolean} , 0) ; 
sector6=xfix({xlBoolean}, 1) ; 
z=xf ix ({xlUnsigned, 10,0}, tempv-5’^ 171) ; 


ramplmod.m 

function z = ramp2(x) 

gain=xfix({xlSigned,20,19} , 1/1800) 

z=xf ix ({xl Signed, 14,13},x’^ gain) ; 


thetaconv2.m 

function [y] = thetaconv(x) 
gainl = xf ix ({xlSigned, 14,10 } , . 14 ) ; 

gain2 = xfix({xlSigned,14,10},1/gainl) 
if x<0 

y=xfix ({xlUnsigned, 10,0}, (x+gainl) ’^gain2’^1024) ; 
else 

y=xf ix ({xlUnsigned, 10,0},x’^gain2’^1024) ; 
end 


C. MATLAB M-FILES USED FOR ENCODER 
mcodeSA.m 

function [NewRotation, NewValue, Direction] = mcode5A(Pk, Nk, Zk, 
OldRotation, OldValue) 

%This faction records CW and CCW rotations via matlab MCODE simulation 
%block. it also records 360 degrees of rotation and total degrees 
traveled 

%Author: LT Andrew M LaValley 
%Last Modified: 14 APR 08 

if Pk && ~Zk && ~Nk %lncrements the # of pulses produced in the CW 
rotation 


Direction=0; %indicates positive direction 

if 01dValue==399 %this keeps the value positive if the value is 
zero. 


NewValue=xfix({xlSigned,15,0},0); 

NewRotation=xfix({xlSigned,15,0},OldRotation+1); 


else 

NewValue=xfix({xlSigned,15,0},OldValue+1); 
NewRotation=xfix({xlSigned,15,0},OldRotation); 
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end 

elseif Nk && ~Zk && ~Pk %Decrements the # of pulses produced in the CCW 
rotation if not zero 

Directional; %indicates negitive direction 

if OldValue==0 %this keeps the value positive if the value is zero. 

NewValue=xfix({xlSigned,15,0}, 399) ; 

NewRotation=xfix({xlSigned,15,0},OldRotation-1); 

%Direction=l; 

else 

NewValue=xfix({xlSigned,15,0},OldValue-1); 

NewRotation=xfix({xlSigned,15,0},OldRotation); 


end 


else %01dValue will equal NewValue and OldRotation will equal NewRotion 
if 

%the above if/else statements do not apply. 

Direction=0; %default direction 
NewValue=xfix({xlSigned,15,0},OldValue); 

NewRotation=xfix({xlSigned,15,0},OldRotation); 

end 


D. MATLAB M-FILES FOR CHIPSCOPE INTEREACE 
Black_box_dc_machine2_config.m 


function code_config(this_block) 

% Revision History: 

g, 

o 

% 18-Dec-2008 (15:15 hours): 

% Original code was machine generated by Xilinx's System 

Generator 

% after parsing H:\Docs\work_files\Docs\classes\EC3130_2009\DC 

% machine lab\black_box_dc_machine2.vhd 

this_block.setTopLevelLanguage( 'VHDL' ); 

this_block.setEntityName( 'code' ) ; 

% System Generator has to assume that your entity has a 

combinational 

% feed through; 

% if it doesn't, then comment out the following line: 
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this_block.tagAsCombinational; 


this_block.addSimulinkInport( 'ind' ); 
this_block.addSimulinkInport( 'ila_clock' ); 
this_block.addSimulinkInport( 'ind2' ); 

this_block.addSimulinkOutport( 'outd' ) ; 
this_block.addSimulinkOutport( 'open_loop' ); 
this_block.addSimulinkOutport( 'speed' ); 

outd_port = this_block.port (' outd ') ; 
outd_port.setType( 'UFix_l_0' ); 

open_loop_port = this_block.port( 'open_loop' ); 
open_loop_port.setType( 'UFix_l_0' ); 
speed_port = this_block.port (' speed ') ; 
speed_port.setType( 'UFix_8_0' ); 

a 

o 

if (this_block.inputTypesKnown) 

% do input type checking^ dynamic output type and 
% this code block. 

if (this_block.port (' ind '). width 1); 

this_block.setError (' Input data type for port 
width=l.' ); 
end 

this_block.port( 'ind' ).useHDLVector(false); 

if (this_block.port (' ila_clock '). width 1); 

this_block.setError (' Input data type for port 
have width=l.' ); 
end 

this_block.port( 'ila_clock' ).useHDLVector(false); 

if (this_block.port (' ind2 '). width 48); 

this_block.setError (' Input data type for port 
width=48. 'Y; 
end 

end % if (inputTypesKnown) 

a 

o 


if (this_block.inputRatesKnown) 

setup_as_single_rate(this_block, 'elk' , 
end % if(inputRatesKnown) 

o, 

o 


ce ' ) 


% Add addtional source files as needed. 
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generic setup in 


"ind" must have 


"ila clock" must 


"ind2" must have 






% I- I Add files in the order in which they should be 

% compiled. | If two files "a.vhd" and "b.vhd" contain the entities 
% entity_a and entity_b, and entity_a contains a | component of type 
% entity_b, the correct sequence of | addFileO calls would be: | 

% this block.addFile('b.vhd'); | this block.addFile('a.vhd'); 


% this_block.addFile(''); this_block.addFile(''); 
this_block.addFile( 'black_box_dc_machine2.vhd' ); 

return; 


function setup_as_single_rate(block,clkname,cename) 
inputRates = block.inputRates; 
uniqueInputRates = unique(inputRates); 

if (length(uniqueInputRates)==1 & uniqueInputRates(1)==Inf) 

block.setError (' The inputs to this block cannot all be constant.'); 
return; 
end 

if (uniqueInputRates(end) == Inf) 
hasConstantInput = true; 

uniqueInputRates = uniqueInputRates(1:end-1); 

end 

if (length(uniqueInputRates) 1) 

block. setError (' The inputs to this block must run at a single 
rate.' ); 

return; 

end 

theInputRate = uniqueInputRates(1); 
for i = 1:block.numSimulinkOutports 

block.outport(i).setRate(theInputRate); 

end 

block.addClkCEPair(clkname,cename,theInputRate); 
return; 
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APPENDIX C: SIMULINK/ XILINX MODEL OF WIND ENERGY 

CONVERSION SYSTEM 


DbJ Fsd Ind S 



H:'ltieEia3VObl Fed Ind 6 mdl 


printed 13-May-2Q09 18 40 


page li'M 
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printed 13-May-2(}09 1S-40 
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APPENDIX D: TRANSFORMATION DERIVATION 


^ab ^bc ~^b~^c~^b~ ’ whcrc = 0 for a wye connected capacitor bank (with the neutral floating) 

Note also that i^ + 4 + 4 “ ^ ^ wye connected transformer winding set. 
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